Alternative model of the Antonov problem 
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Abstract 

Astrophysical systems will never be in a real Thermodynamic equilibrium: they undergo an evaporation process due to the 
fact that the gravity is not able to confine the particles. Ordinarily, this difficulty is overcome by enclosing the system in a 
rigid container which avoids the evaporation. We proposed an energetic prescription which is able to confine the particles, 
leading in this way to an alternative version of the Antonov isothermal model which unifies the well-known isothermal and 
polytropic profiles. Besides of the main features of the isothermal sphere model: the existence of the gravitational collapse 
and the energetic region with a negative specific heat, this alternative model has the advantage that the system size naturally 
appears as a consequence of the particles evaporation. 

PACS numbers: 05.20.-y, 05.70.-a 
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I. INTRODUCTION 

Thermodynamical properties of selfgravitating systems 
are very different from the ones exhibited by the tradi- 
tional systems. They are typical nonextensive systems 
since they are nonhomogeneous and the total energy is 
not extensive, which is a consequence of the long-range 
character of gravitational interaction. They also exhibit 
energetic regions with a negative specific /leai, whichper- 
sist even in the thermodynamic limit [ISIEIlHIElEll- 
That is the reason why the Gibbs canonical ensemble is 
non applicable to the description of selfgravitating sys- 
tems, since this ensemble is not able to access to those 
macroscopic states possesing a negative heat capacity. 
At first glance, the selfgravitating systems could only be 
described by using the microcanonical ensemble. 

Gravitation is not able to confine the particles: it is 
always possible that some of them have the sufficient en- 
ergy for escaping out from the system, so that, the self- 
gravitating systems always undergo an evaporation pro- 
cess. Therefore, they will never be in a real thermody- 
namic equilibrium. This difficulty is usucdly avoided by 
enclosing the system in a rigid container |a, la, 13 1 which 
could be justified when the evaporation rate is small and 
certain kind of quasistationary state might be reached. 
There are some approaches which have taken into ac- 
count this kind of regularization of the long-range sin- 
gularity of the Newtonian potential in a microcanonical 
framework |^[To|. 

The use of the rigid container can be conveniently sub- 
stituted by imposing the following energetic prescription: 
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-Hm(/)(r) < es < 0, 



(1) 



where ^P^ is the kinetic energy of a particle of mass m 
and mcj) (r) its gravitational potential energy at the point 
r, being es an energy cutoff which is determined from the 
existence of certain tidal forces. All those particles sat- 
isfying the condition ^ are confined by the gravity, on 
contrary, they will be able to scape out from the system 
if they do not lose their excessive energy. This regulariza- 
tion procedure is characteristic of the Michie-King model 
for globular clusters 11], which supposes that those stars 
that gain suficient velocity through encounters are able 
to escape or are removed by tidal forces. The main mo- 
tivation of such regularization scheme relies on the con- 
sideration of evaporation effects, and therefore, this reg- 
ularization procedure is more realistic than the box reg- 
ularization (the use of the rigid container). The aim of 
this paper is to develop an alternative model to the stan- 
dard isothermal sphere model of Antonov 1^ based on 
the consideration of this energetic prescription starting 
from microcanonical basis. 



II. STATISTICAL DESCRIPTION 

Let us consider the iV-body selfgravitating Hamilto- 
nian system: 
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Taking into consideration the above regularization pre- 
scription the admissible stages of this system are 
those in which the kinetic energy of a given particle sat- 
isfies the condition: 
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pl<Uk = m[0s - 0(rfc)] , 



(3) 
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being (/) (r^) the gravitational potential at the point Vk 
where the fc-th particle is located: 
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(4) 



where we introduce the tidal potential 0s (eg — m(j)s). 

The regularized microcanonical accessible volume Wr 
is given by: 



Wr^^ J^ S[E-HN]dX, 

where Xji is a subspace of the iV-body phase-space 
where the condition J^Jl takes place, being dX = 

^3W^^3Arp/(27rn)3^ = II^^tSS^ t^ie volume element. 

k 

The spacial coordinates should be also regularized in or- 
der to avoid the short-range divergenge of the Newtonian 
potential, which will be performed below by using the 
mean field approximation. This integral is rewritten by 
using the Fourier representation of the delta function as 
follows: 
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where Zf{ (z, A^): 

ZRiz,N)^l^ J^ e^p{-zHN)dX, (6) 

is the canonical partition function with complex argu- 
ment z — f3 + ik , with /3 G 5R. Integration by d^^P 
yields: 



If ,3Ar„/ TO 



NlJ^3N \2TTh?Z 

being x(i?;z): 



exp [-zUn [R) + X (i?; 2:)] , 

(7) 
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:(i?;z) = ;^lnF[yi^. 



(8) 
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where F (z) is defined by: 



F (z) = erf (z) -=z exp (— z^) 



(9) 



and shown in the figure The asymptotic behaviors 
of F (z) are given by: 

^, ^ r ^^'■^ + O (^^) when z ^ 

F(z)^ 37. ^^^^^ . (10) 
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FIG. 1: Representation of the function F [z). Its symptotic 
behaviors are also shown. 



We are interested in describing the large N limit. This 
aim can be carried out by using following procedure: we 
partition the physical space in cells {cq.} being Va the 
positions of their centers. We denoted by Ua — n{ra) 
the number of particles inside the cell at the position . 
Introducing the density p (r^) = n {va) /va, being Va the 
volume of the cell c^, the functions {R) and x {R'l 
are rewritten by using a mean field approximation as fol- 
lows: 
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(11) 



X (R; z) — > X [p, 0; z] 

= J^^ p (r) InF (v/zm [cps - 4> (r)]) d^r, (12) 
being (r) the Newtonian potential for a given p profile: 



^{v)^g[p] = - j 



Gmp {ri) d^ri 



/g?3 |r-ri| 
which is the Green solution of the problem: 

= 4nGp. 



(13) 



(14) 



It is easy to understand that this mean field approxi- 
mation acts as a partial regularization procedure for the 
short-range singularity of the Newtonian potential. The 
microscopic fluctuations of the Newtonian potential are 
disregarded in this approximation, since this field is con- 
sidered constant inside the volume of each cell. Thus, the 
gravity effects on the microscopic picture of each cell are 
reduced to the truncation of the velocity distribution of 
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the particles. This regularization is partial since it does 
not avoid the gravitational collapse 

Using the partition of the physical space in cells, the 
integration by d'^^i? can be approximately given by: 



where 



(15) 



Vp (r) V<i> (r) b {4> (r) -Q{p\}b{N ~N {p\) x 



X exp {-p/ [p; z\ - zU [p, (f>] + x [p, <!>■, z]} ■ (23) 

The delta functions are also conveniently rewritten by 
using their Fourier representation: 



(^) = i,tL:wise -dE-EE--- (16) 



The following factor can be rephrased in the mean field 
approximation as follows: 
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being 



Pf [P: - 



lnp(r)-l + -ln 



(18) 

where the Stirling formula Inn! ~ nlnn — n was used. 
Finally, we rephrase the summation in the mean field 
approximation as follows: 



Taking into consideration all approximations intro- 
duced above, the canonical partition function @ is 
rewritten as: 



Zc[z,N] = Jvp{r)5[N-N[p]] x 

X exp {-Pf [p, z] - zU [p, (j> [p]] +x[p,(f> [p] ; z]} , (20) 
being N [p] the particles number functional: 



N[p]^ [ d^rpir) 



(21) 



In order to avoid the complicate p dependence of ((> [p] in 
(|20|l . we introduce the identity: 



(22) 



V^{r)S{cf>{r)-g[p]}^l, 
in the functional integral (|20|l : 



S {cj, (r) - g [p]} ^ / m (r) exp [ / d'rJ (r) {<j> (r) - ^ [p]} 

(24) 



as well as 
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6{N-N[p])^ /■^exp[zi(iV-7V 
J 27r 



(25) 



being zi — p + iq with /i G SR; J (r) = j (r) + ih (r), a 
complex function with j (r) and h (r) £ 5R. Thus, the 
canonical partition function Zc [z, N] is finally expressed 
as: 
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/ S /2'P(r)I^<^(r)I^/i(r)exp{ziiV-H[p,0;z,zi,J]}. 

— OO 

(26) 

The functional H [p, (jj; z, zi, J] is given by: 



= -P/ [p; ^] + [P> (t>]-X [p, z] + ziN [p] + if [0, p, J] , 

(27) 

being K [0, p, J] the exponential argument of the expres- 
sion if^ . 

The reader may check that when N is scaled as ^ 
aN , and the following quantities are scaled as follows: 



p a p, V ^ a 3r, 0^Q!3 0, 
z — > a 3z, zi— >zi, J-^aiJ, 



(28) 



all terms in 127|l scale proportional to a, and therefore: 



H [p, 0; z, zi, J] ai? [p, 0; z, Zi, J] . 



(29) 



The thermodynamic limit is carry out tending a to the 
infinity, a — > oo. Thus, we can estimate Zc [z, N] for N 
large by using the steepest decent method. The Planck 
potential V (J3, N) = — In Zc [(3, N] is thus obtained as 
follows: 



V [(3, N] ~ -max min [pN - H [p, 0; p, j]] . (30) 
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where the stationary conditions: 
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(31) 



lead to the foUowing relations: 
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-/i - -/3m0 + g [j] 



j = - -j^Pmp + pdtj, In F ^ jim [c^s 



= g [p] and TV = iV [p] , 



(32) 



(33) 



(34) 



The relations (|32|l and H33|) define the state equation; 
and the relations (|34|l establish the Newtonian potential 
for a given p profile as well as the normalization con- 
strain for the number of particles. The relation H32|) is 
conveniently rewritten by using (|33|l as follows: 

p = iV(/3,M)exp[$ + C]i^($5^ , (35) 

where $ = /3m \(^s ~ <f\ and N {f3, p) — 

(2;;^) '^^P ~ ^"^'^•s), being C a new func- 
tion which is given by: 



C = -I3mg pa* In F $2 



(36) 



It is not difficult to show that the Planck potential is 
given by: 



(37) 

The Boltzmann entropy can be estimated by using the 
steepest decent method as follows: 



SBiE,N)c^mm[PE-V[P,N]], (38) 
being E [(3, p, cj)] the energy functional: 



3 -$9$ In ($3^ 



It can be easily seen that the quantity: 
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e(/3,$) = 3 - $5<i,lnF ($5^ 
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(40) 



represents the kinetic energy per particle at a given point 
of the selfgravitating system. Note that p vanishes when 
$ tends to zero in lp!K|l . so that, the particle distribution 
has been regularized. This state equation differs from 
the one obtained in the isothermal model by the pres- 
ence of the truncation function as well as the 



function C. Although C naturally appears in 
our derivation, its existence is closely related with the 
modification provoked in the microscopic picture of the 
system by the evaporation. An example of this affirma- 
tion is the deviation of the Maxwell distribution along 
the system, which can be noted in the kinetic energy per 
particle e (/?, <I>) (|40|l . A naive energy truncation of the 
Maxwell- Boltzmann distribution: 



^MB = Co exp 
leads to the state function: 



13 1 — + mc, 
2m 



(41) 



P' 



exp [4>] ^$2^ 



(42) 



but here the driving function C does not appear. 

According to the scaling laws H28|l . this astrophysical 
model obeys to the following thermodynamic limit: 



N 00, keeping constant — - and LN^ , (43) 

being L the characteristic system size. This thermo- 
dynamic limit was obtained in ref.(12| for the self- 
gravitating fermions model. The reader may surprise of 
the iV-dependence of the characteristic system size L. 
However, there is nothing strange in this behavior since 
the selfgravitating gas is constituted by punctual parti- 
cles, and therefore, this model system can be reduced to 
a point in the thermodynamic limit — > 00. 

It is easy to understand that this is an asymptotical be- 
havior which disappears when the particles size or/and 
the relativistic effects are taken into account. It means 
that it should exist a superior limit of N in which these 
A-dependences of the energy and system size become 
invalid, and therefore, the selfgravitating nonrelativistic 
gas model is inapplicable for describing the thermody- 
namical properties of such massive astrophysical systems. 

It is interesting to analyze the asymptotic behavior of 
the particles density state function p (<I>; /3, p). According 
to the asymptotic behaviors H1U|I of the function F (z), 
the state equation (|35|l becomes in the isothermal distri- 
bution when $ > 2.5 : 



p oc exp [$] 



(44) 



which is characteristic of the inner regions of the system 
(the core). At local level, the particles velocities obey to a 



4 



Maxwell distribution, where the kinetic energy per parti- 
cle H4()ll is given by 3/2/3. When we move from the inner 
regions towards the outer ones, the kinetic energy per 
particle decreases until zero, being this approximately 
given by e ~ 3<i>/5/3. There the particles velocities obey 
to an uniform distribution, and so, the state equation in 
the halo obey to a polytropic distribution: 



p oc 



30F 



$ 2 



(45) 



All these behaviors appear as consequence of the high 
energy cutoff (PJ in the Maxwell distribution for the ki- 
netic energy / (fc) = 27r~2 fcs exp (— fc) {k = /3p^/2m ), 
which is shown in the figure (|2Jl- Note that the cutoff 
value of k is equal to In the inner regions where 
<f> > 2.5, this energy cutoff is unimportant because of 
it is located in the tail of the distribution. However, 
it modifies considerably the character of the state func- 
tion for p in the halo due to the fact that it is located 
before the pick. The particles density p will obey to the 
polytropic dependence (|45(l throughout the whole system 
volume when $ is also small in the inner regions. 



f(k) 




FIG. 2: Maxwell distribution for the kinetic energy. The char- 
acter of the state function for the particles density p depends 
on the localization of the energy cuttof "I>. 



Thus, the consideration of the energetic prescription 
always leads to the existence of a polytropic halo, 
where the system core can be isothermal or polytropic, 
according to the values of the function $ in the inner re- 
gions of the system. Due to the general properties of the 
spherical solutions with a polytropic profile, the particles 
density will vanish at a finite radio Rudaii which can be 
identify with the characteristic size of the system. This 
radio is related with the tidal potential (ps throughout 
the relation: 



being M the system total mass. Thus, the characteristic 
size of the system is determined by the tidal interactions. 



III. NUMERICAL STUDY 

In order to perform a numerical study of the equations 
obtained in the previous section, we express the energy 
in units of Eq — GA'P / Rudai and the lenght in units of 
the system size Rudah and the mass in unit of M. It is 
not difficult to show that the functions $ and C obey to 
the following structure equations: 



Ar<i> = -AnFi (^', $) , Ar"^ = -AttF2 $) , (47) 

where A^u = r~^dr (r^dru) is the radial part of the 
Laplace operator, being the functions Fi (5*, $) and 
F2 (5', $) defined by: 



Fi (*, $) = exp {^ + ^)F (^$5^ 



F2 (*, $) = —;= exp (^) $2 . 
a/tt 



(48) 



In the expressions above the function '5 = C+ln/C (/3, /i), 
with K, {/3, fjL) = /3~2 exp (— /j, -I- /?). The solutions of the 
nonlinear system H47(l satisfy the following conditions at 
the surface: 



$ (1) = 0, (1) = C (1) = -C" (1) 



(49) 



which are derived from the boundary conditions of Green 
solutions ifT^ . 

The solution can be obtained from the imposition of 
the following boundary conditions at the origin: 



$(0) = $0 > 0, 
*(0) = ^'($o), 



(50) 
(51) 



where the value of 4* (0) depends on the parameter $0 
because of $ must vanish when r = 1 . This situation can 
be overcome redefining the problem as follows: 



l'(r) = V(e)+21ni?,„, $(r) =^(0 



(52) 



being {tp (^) , ip (^)} the solution of (|47|) whose boundary 
conditions are given by: 



(0) = $0 > 0, if' (0) = 0, ^ (0) ^ 0, ij' (0) = 0, (53) 
where f is related with r throughout the relation: 



GM 

Rtidal ' 



(46) 



(54) 
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being Rm the radio at which (p vanishes {(p {Rm) = 0). 

The canonical parameter f3 and are obtained from 
the relations: 

i?„aj(^(i?„0 = -/3 = -/3(*o), (55) 



+ --/i, (56) 

being h — h ((f>o) = — In [/C (/?, /i) /i?^„] , which allow us 
to express as a function of $0^ 

/i = /i + /3- iln/3-21ni?,„. (57) 

Taking into consideration the equations (|39|1 and (|37|1 . 
the total energy and the Planck potential per particle are 
rewritten as follows: 

V (<I>o) = -1 - A* + i/? + /i ($0) + -^h2 ($0) . (59) 
being 

/ii ($0) = 4^e'rfe [^1 ^) + ^2 (^, <p)] , (60) 

($0) - J^"""^ ^T^ed^ (v., p) + . (61) 

IV. RESULTS AND DISCUSSIONS 

Figures JSJ and Q show respectively the caloric curve 
and the central density of the model. These dependences 
evidence that this model system exhibits the main fea- 
tures of the Antonov isothermal model: the existence of a 
negative specific heat when < e < , and the gravita- 
tional collapse for e < eA, where the central density grows 
towards the infinity and the system develops a core-halo 
structure, being ca = —0.806 and e_B = —0.446. 

This conclusion is supported by the analysis of the 
thermodynamical potentials of the model: the entropy 
in the microcanonical ensemble, at figure ((SJ, and the 
Planck potential in the canonical ensemble, at figure 
©. The figure (|SJ) shows that the points of the superior 
branch of the caloric curve ||2J) correspond to equilibrium 
configuration while the others represent unstable saddle 
points. No equilibrium states exist when e < eA- 

On the other hand, the figure © evidences that the 
canonical ensemble can access only to those equilibrium 
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FIG. 3: The alternative model exhibits the main features of 
the isothermal model of Antonov. 



1 E- 

1 r 



10^ r 




10° r 



-101 2 3 4 5 

energy 

FIG. 4: Central density versus energy: this dependence shows 
the existence of the gravitational collapse for e < eA- 



states beloging to the interval < (3 < (3b (with e > eb), 
where /3b = 1-193. The energetic region < e < e_B 
is invisible for the canonical description due to the neg- 
ativity of the heat capacity. No equilibrium states exist 
for 13 > (3b- This fact evidences the existence of a grav- 
itational collapse in the canonical ensemble beyond the 
critical point 0b , which is usually refered as an isother- 
mal collapse [13 • 

The Gibbs' argument, the equilibrium of a subsystem 
with a thermal bath, is non applicable to this situation 
because of no reasonable thermal bath exits for the astro- 
physical systems. Therefore, the isothermal catastrophe is 
not a phenomenon with physical relevance since it can be 
never obtained in nature: the consideration of a thermal 
bath in the astrophysical system is outside of context. A 
different significance possesses the gravothermal catas- 
trophe. The gravitational collapse is the main engine of 
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with a linear dimension R = Rudai 



FIG. 5: Entropy versus energy: Equilibrium configurations 
exit only for e > ea. 
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FIG. 7: The caloric curves. The figure shows the comparison 
between the two models. 
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FIG. 6: Planck potential versus inverse temperature: the 
canonical ensemble access only to those equilibrium stages 
with /? < /3b (e > es). 



structuration in astrophysics and it concerns almost all 
scales of the universe: the formation of planetesimals in 
the solar nebula, the formation of stars, the fractal nature 
of the interstellar medium, the evolution of globular clus- 
ters and galaxies and the formation of galactic clusters 
in cosmology [isj . 

Let us now carry out a comparative study between the 
isothermal model and this alternative onel. As already 
mentioned, the only difference between these approaches 
relies on the regularization prescription of the long-range 
singularity of the selfgravitating gas: the isothermal 
model avoids the particles evaporation by using a rigid 
container, while the present model takes into account the 
effect of this evaporation by truncating the kinetic energy 
distribution function. This study is carried out by con- 
sidering a spherical container in the isothermal model 
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FIG. 8: The central density versus energy. The central den- 
sity is much greater in this alternative version of the Antonov 
problem. 

Figure J?)) shows the caloric curves of these models. In 
spite of the qualitative similarity of these dependences, 
the alternative model is able to describe an additional en- 
ergetic range: from —0.806 to —0.335; but the isothermal 
model describes equilibrium configurations belonging to 
the interval 1.193 < (3 < 2.518, which are cooler than the 
ones described by using the first model.. 

A second difference is evidenced in regard to the cen- 
tral density versus energy dependence, which is shown 
in the figure JH)). The central density at the critical en- 
ergy of the gravitational collapse is much greater in the 
alternative model than the isothermal one, and this qual- 
itative relationship seems to be applicable to the whole 
energetic range. 
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isothermal radio 
mass content 




FIG. 9: Isothermal core radio and its mass content in the 
alternative model. 



The figure ^ shows two interesting observables: the 
radio in which the system exhibits an isothermic behav- 
ior, Ri,c, and the mass content enclosed inside, M^ c- 
This quantities characterizes the size and the mass of 
the isothermal core. The reader may observe that the 
isothermal core contains at the critical energy of the grav- 
itational collapse, the half of the system total mass inside 
the ~ 1% of the system volume. On the other hand, it 
is interesting to note that this isothermal core disappears 
at e* ~ -0.066. 

The quantitative differences in the thermodynamical 
description between these models seem to be explained 
by the following reasonings. Most of the energy con- 
tribution to the total energy comes from the isothermal 
core. The contribution of the gravitational potential en- 
ergy is dominant in the core due to the high mass concen- 
tration enclosed inside this region, where moreover, the 
Newtonian potential exhibits its highest values. Outside 
the isothermal core, the kinectic energy contribution de- 
creases considerably as consequence of the deviation from 
the isothermal character of the microscopic particles dis- 
tribution function due to the evaporation. 

These arguments can be rephrased as follows: in the 
isothermal region, where —0.806 < e < —0.066, the 
present model behaves as an isothermal model with a 
characteristic size equal to the size of the isothermal core. 
Since the relevant canonical variables in the isothermal 

model is ?7 = PGM/R mill, a rough estimation of the 

maximal value /3max — Pb in which the isothermal col- 
lapse takes places is given by: 



R 



M; 



t-c nisoth. 
max 



(62) 



By using the characteristics values Ri,c ~ 0.2, M^.c ~ 0.5, 
and /3max''' ~ 2.5 (obtained from the isothermal model), 
the formula H62|) gives a fairly good estimation of /3b w 1 
{J3b = 1.193). 



The main difference between these model is in regard 
to the character of the equilibrium profiles in the high en- 
ergy region with e > 0. As already discussed, the isother- 
mal core has disappear in this energetic region and the 
equilibrium configurations of the system are essentially 
polytropic, while the isothermal model leads to a uni- 
form distribution of the particles throughout the volume 
of the rigid container. 

In this case, the structure equations H47() become in a 
quasi-polytropic model: 



A$ = -47r exp (*) — 



A* 



-47rexp (*) 



(63) 



which exhibits the same fractal characteristics of a poly- 
tropic model with polytropic index 7 = f • This poly- 
tropic index characterizes an adiabatic process of the 
ideal gas of particles, which is in our case the evapo- 
ration of the system in the vacuum. However, this equa- 
tion system is not equivalent to the polytropic model due 
to the presence of the driving function C. This purely 
polytropic profile is obtained by disregarding the second 
equation and setting 5" = in the first one. 




FIG. 10: Comparison among the equilibrium profiles; A) Al- 
ternative model with an isothermal core, B) Isothermal pro- 
file, C) Alternative model with a quasi-polytropic profile, and 
D) Purely polytropic profile. 



Let us finalize our discussion by carrying out a com- 
parative study among the equilibrium profiles obtained 
from the alternative model with the well-known isother- 
mal and polytropic profiles. These equilibrium profiles 
are shown in the figure (jlO|l . Equilibrium profiles A and 
C were obtained from the alternative model: profile A 
corresponds to an equilibrium configuration possessing 
an isothermal core, while C is a quasi-polytropic equilib- 
rium profile. Profiles B and D correspond respectively to 
an isothermal and polytropic configurations. The reader 
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can note that the alternative model profiles differ essen- 
tially in regard to the existence of the isothermal core, 
since no qualitative differences are evidenced in the halo 
structure. On the other hand, these results suggest that 
a system undergoing an evaporation process concentrates 
their particles in the inner regions more than what the 
isothermal or the polytropic models predict. 

V. CONCLUSIONS 

As already shown in the present paper, the considera- 
tion of the energetic prescription leads to an alterna- 
tive version of the Antonov problem. This model, besides 
of exhibiting the main features of the isothermal model: 
the gravitational collapse and the energetic region with 
a negative specific heat, has the advantage that the sys- 
tem size naturally appears as consequence of the particles 
evaporation. There is no need of enclosing the system in 
a rigid container in order to avoid long-range singularity 
of the Newtonian potential because this regularization 
procedure is sufhcient to access to finite equilibrium con- 
figurations. 

It is remarkable that the present approach unifies the 
well-known isothermal and polytropic equilibrium pro- 



files. As already mentioned, the equilibrium profiles de- 
rived from this alternative model differ essentially in re- 
gard to the existence of the isothermal core, since no 
qualitative differences are evidenced in the halo struc- 
ture. The comparative study of the equilibrium p profiles 
suggests that a system undergoing an evaporation pro- 
cess concentrates the particles in the inner regions more 
than what predict by the isothermal or the polytropic 
models. 

There are many open questions in the study of this 
model system by using this kind of regularization scheme, 
as example, the dynamical aspects: Which is the infiu- 
ence of the particles evaporation in the system dynam- 
ical evolution? and how could it be performed this dy- 
namical description by preserving this energetic prescrip- 
tion? Further studies should clarify these questions. The 
present study is carried out by considering that the evap- 
oration rate is small enough in order to ensure that the 
system reaches a quasistationary state. This assump- 
tion is satisfied by the globular clusters, since the relax- 
ation time treiax in them differs considerable from the 
evapotation time tevap, treiaz ^ tevap (see in ref._14]). 
The dynamical approch could be developed by consid- 
ering that the system evolutes slowly throughout these 
quasistationary states. 
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